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PREFACE 

The  work  reported  herein  was  conducted  by  the  Arnold  Engineering 
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Program  Element  65807F.  The  results  were  obtained  by  ARO,  Inc.,  AEDC 
Division  (a  Sverdrup  Corporation  Company) , operating  contractor  for  the 
AEDC,  AFSC,  Arnold  Air  Force  Station,  Tennessee,  under  ARO  Project  No. 
P33A-G4A.  The  author  of  this  report  is  David  L.  Whitfield,  ARO,  Inc. 
Elton  R.  Thompson  is  the  Air  Force  project  manager.  The  manuscript 
(ARO  Control  No.  ARO-PWT-TR-77-43)  was  submitted  for  publication  on 
June  16,  1977. 

The  author  would  like  to  acknowledge  Dr.  J.  L.  Jacocks,  ARO,  Inc., 
with  whom  periodic  discussions  were  held  during  the  course  of  this  work 
which  contributed  appreciably  to  improving  the  end  result. 
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1.0  INTRODUCTION 

Thi9  report  is  concerned  with  the  development  of  an  analytical 
expression  for  the  velocity  distribution  in  a turbulent  boundary  layer. 
Because  numerous  velocity  distribution  expressions  have  been  developed 
for  various  portions  of  a turbulent  boundary  layer  (see,  e.g..  Ref.  1), 
a summary  of  the  more  recently  developed  or  frequently  used  methods 
follows  in  order  to  place  the  present  investigation  in  proper,  perspective. 

A commonly  used  expression  for  the  turbulent  boundary- layer  veloc- 
ity distribution  is  the  logarithmic  distribution  or  law-of-the-wall. 

The  logarithmic  velocity  distribution  is  valid  in  a region  near  the  wall 
but  not  in  the  immediate  vicinity  of  the  wall  in  the  so-called  viscous 
sublayer  and  buffer  layer,  nor  is  it  valid  in  the  outer  part  of  the 
boundary  layer  in  the  so-called  wake  region.  Particular  expressions  for 
the  sublayer  and  buffer  layer  have  been  developed  (Ref.  1),  and  a single 
expression  for  the  sublayer,  buffer  layer,  and  logarithmic  region  was 
developed  by  Spalding  (Ref.  2).  For  the  outer  portion  of  a turbulent 
boundary  layer.  Coles  (Ref.  3)  developed  the  law-of-the-wake.  Since  the 
appearance  of  Coles'  work,  other  expressions  for  the  law-of-the-wake 
were  developed  which  correct  certain  deficiencies  of  Coles'  (Ref.  3) 
law-of-the-wake  as  pointed  out  by  Dean  (Ref.  4).  Both  Spalding’s 
expression  (Ref.  2)  and  the  law-of-the-wake  were  developed  such  that 
they  asymptote  to  the  logarithmic  velocity  distribution.  Dean  (Ref.  4) 
used  this  feature  of  commonality  in  asymptotic  behavior  to  combine 
Spalding's  (Ref.  2)  expression  and  a law-of-the-wake  to  obtain  one 
unified  expression  that  is  valid  from  the  wall  to  the  outer  edge  of  the 
boundary  layer.  Dean's  expression,  however,  cannot  be  written  such  that 
either  the  velocity,  u,  or  the  distance  from  the  wall,  y,  is  an  explicit 
function  of  the  other.  This  is  a consequence  of  using  Spalding's  (Ref. 

2)  expression  which  gives  y+  explicitly  in  terms  of  u+.  The  law-of-the- 
wake  used  by  Dean  (Ref.  4)  contains  Coles'  profile  parameter,  FT.  In 
general,  II  is  determined  by  locating  or  choosing  6 (which  Clauser  (Ref.  5) 
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calls  an  Ill-defined  quantity)  and  then  using  the  law-of-the-wake  to 
solve  for  H;  or,  II  and  6 are  solved  for  simultaneously  using  the  law-of- 
the-wake.  For  example,  in  the  matching  of  Coles'  law-of-the-wake  as 
given  in  Ref.  6,  with  the  data  presented  in  the  Stanford  Conference 
Proceedings  (Ref.  6),  both  IT  and  <5  were  used  as  free  parameters.  The 
values  of  5 given  in  Ref.  6 varied  considerably  from  values  corresponding 
to  the  usual  definition  of  6,  where,  for  example,  6 is  the  value  of  y 
where  u/ug  = 0.99.  Therefore,  based  on  this  brief  summary  of  available 
analytical  techniques,  a contribution  in  this  area  would  be  to  analyti- 
cally describe  turbulent  boundary- layer  velocity  distributions  by  an 
expression  that  gives  the  velocity  explicitly  as  a function  of  distance 
from  the  wall,  and  that  contains  parameters  which  are  explicitly  defined 
and  can  be  calculated  in  terms  of  properties  of  the  boundary  layer.  The 
purpose  of  this  report  is  to  develop  such  an  expression  and  validate  it 
by  comparisons  with  experimental  data.  The  resulting  expression  should 
be  useful  for  making  hand  calculations  of  complete  turbulent  boundary- 
layer  velocity  profiles,  providing  initial  conditions  for  numerical 
solutions  of  the  boundary- layer  equations,  and  for  use  in  integral 
boundary-layer  theories. 

Attention  here  is  focused  on  two-dimensional  turbulent  boundary 
layers  on  smooth  impermeable  walls.  The  approach  taken  is  the  con- 
struction of  a composite  expression  in  terms  of  the  inner  variable  y+ 
and  the  outer  variable  y/9.  An  analytical  solution  of  the  turbulent 
boundary-layer  equations  (including  a turbulent  kinetic  energy  equation) 
for  the  region  near  the  wall  (Ref.  7)  is  used  for  the  inner  solution, 
and  an  empirical  expression  is  derived  for  the  outer  region  of  the 
boundary  layer.  The  analytical  development  is  given  In  the  following 
section  and  the  accuracy  of  the  analytical  result  is  investigated  in 
Section  3.0  by  comparisons  with  a variety  of  experimental  data,  which 
include  boundary  layers  near  separation,  a boundary  layer  downstream  of 
reattachment,  nonequilibrium  boundary  layers,  and  compressible  boundary 
layers.  A summary  of  the  procedure  and  equations  for  the  computation  of 
turbulent  boundary-layer  velocity  distributions  is  given  in  Table  1. 
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2.0  ANALYTICAL  DEVELOPMENT 


The  derivation  of  the  analytical  result  is  contained  in  this  section. 
The  development  is  carried  out  with  all  variables  in  incompressible 
form.  In  the  following  section,  comparisons  are  made  with  compressible 
data  and  at  that  point  the  variables  will  be  transformed  to  account  for 
compressibility. 


The  desired  composite  expression  is  required  to  have  the  following 

properties:  (1)  recover  the  solution  developed  in  Ref.  7 for  the  inner 

+ 2 + 
region,  i.e.,  y < 0(10  ),  (2)  approach  the  proper  limiting  value  of  u 

-*■  u*  (or  u/ufi  -*■  1)  as  y -*■  »,  and  (3)  recover  the  velocity  profiles 
similar  to  those  correlated  by  von  Doenhoff  and  Tetervin  (Ref.  8)  away 
from  the  wall  in  the  outer  variable  y/9.  Requirement  number  (1)  com- 
prises the  inner  solution  denoted  as  and  requirements  (2)  and  (3) 

+ + 

comprise  the  outer  solution  denoted  as  uq . The  composite  solution  u 
is  taken  as 


u+  = 


“i  + 


(1) 


Consideration  is  now  given  to  particular  expressions  for  u*  and  u*. 


The  derivation  of  the  inner  solution  is  given  in  detail  in  Ref,  7. 
It  was  obtained  from  an  analytical  solution  of  the  turbulent  boundary- 
layer  equations  including  a turbulent  kinetic  energy  equation  which  was 
developed  in  Ref.  7 for  turbulence  modeling.  The  primary  assumptions 
made  in  obtaining  the  inner  solution  were: 


1.  the  total  shear  stress  consisted  of  the  molecular  and  turbu- 
lent contributions,  was  constant,  and  equal  to  the  wall  value. 


2.  the  density  was  constant  and  equal  to  the  wall  value,  p^. 
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3.  the  turbulent  kinetic  energy  was  proportional  to  -u'v',  and 


4.  the  turbulent  kinetic  energy,  or  -u’v',  was  an  explicit  func- 
tion of  u . 


Assumption  (1)  is  commonly  used  in  the  region  near  the  wall  and  assump- 


tion (2)  was  used  to  satisfy  r = (u3u/3y) 

w y 


because  u+  = y+  when  y -> 


0.  Assumption  (3)  is  not  an  uncommon  assumption  in  the  fully  turbulent 
part  of  a boundary  layer  (see,  e.g..  Ref.  7 for  some  history  of  the  use 
and  development  of  this  relation);  however,  some  supporting  evidence  is 
desirable  to  justify  this  approximation  for  use  in  the  inner  region. 
Such  evidence  is  given,  for  example,  by  Fig.  5 of  Hanjalic  and  Launder 
(Ref.  9)  where  they  compared  their  numerical  calculations  of  turbulent 


kinetic  energy  and  shear  stress  (-u'v1)  with  experimental  data  in  the 

es 
2 


+ 2 

inner  region  (y  <0(10)).  If  one  traces  their  calculations  of 


2 + 

-u'v'/u^.  as  a function  of  y and  overlays  it  on  their  calculations  of 

turbulent  kinetic  energy,  one  finds  that  the  distributions  are  similar 

except  for  y+  less  than  about  3.  For  y+  < 3,  their  calculations  of 
2 

-utv7ut  fall  below  the  experimental  data.  In  considering  some  of  the 
data  presented  by  Hanjalic  and  Launder  in  Fig.  5 of  Ref.  9,  Hlnze  (Ref. 
10)  observed  that,  whereas  it  can  be  proved  that  -u'v1  must  behave  as 
y+  or  y+  as  y+  + 0,  experimental  data  indicate  that  -u'v*  ~ y+  for 
y as  small  as  measurements  are  currently  available  which  is  about  y+  of 


0(10  ).  This  experimentally  observed  quadratic  variation  of  -u*v’/u 
with  y near  the  wall  is  consistent  with  the  experimentally  observed  and 
numerically  calculated  distributions  of  turbulent  kinetic  energy  near 
the  wall  (see,  e.g..  Ref.  9,  Fig.  5 again).  Therefore,  based  on  these 
experimental  and  numerical  distributions  of  turbulent  kinetic  energy  and 
shear  stress,  assumption  (3)  appears  reasonable  in  the  Inner  region  as 
well  as  the  outer  portion  of  the  boundary  layer  although  the  constants 
of  proportionality  may  differ.  Assumption  (4)  Is  a consequence  of  the 
fact  that  both  u+  and  -u’v*  are  functions  of  y+  in  the  inner  region,  and 
hence  -u’v*  is  a function  of  u+. 
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The  solution  procedure  in  Ref.  7 was  to  solve  the  turbulent  kinetic 
+ 

energy  equation  for  -u'v'  as  a function  of  u (recall  that  the  turbulent 
kinetic  energy  was  assumed  proportional  to  -u'v').  Then  -u'v'  was  used 
in  the  expression  for  the  total  shear  stress  (which  was  assumed  constant 
and  equal  to  the  wall  value)  to  obtain  a first-order  differential 
equation  for  the  velocity.  This  fir9t-order  equation  was  obtained  by 

normalizing  the  expression  for  the  total  shear  stress,  t = p(3u/9y)  - 

+ + 2W 

pu'v’,  which  resulted  in  the  equation  du  /dy  - pu'v’/u  = 1.  This 

1 +T  + 

equation  was  then  solved  with  the  boundary  condition  u.  = 0 at  y =0 

+ + 

to  obtain  the  solution  for  u^  as  a function  of  y (see  Ref.  7), 

The  significance  of  the*  results  of  Ref.  7 is  that  a consistent 
description  of  mean  turbulence  quantities  in  the  inner  region  was  ob- 
tained in  exceedingly  simple  mathematical  form.  This  is  indicated  by 
Figs.  1,  2,  and  3 and  the  corresponding  inner  solutions  for  the  veloc- 
ity, Reynolds  stress,  turbulence  production,  and  direct  dissipation  of 
mean-flow  energy  given,  respectively,  as 

ut  = (0.09v+) 

1 0.09  ' 


(3) 


(4) 


(5) 
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Note  from  Eq.  (2)  that  the  classical  sublayer  result  of  u+  =*  y^"  is 

+ 

recovered  for  small  values  of  y . Also  of  interest  is  the  accuracy  in 
the  slope  of  the  velocity  distribution  as  Indicated  by  the  comparison  of 
the  dissipation  term  [Eq.  (5)J  with  experimental  data  in  Fig.  3.  The 
Reynolds  stress  and  turbulence  production  comparisons  are  included  to 
complete  the  turbulence  modeling  description  according  to  the  inner 
solution  and  illustrate  its  simplicity  and  good  agreement  with  experi- 
mental data. 


In  contrast  to  u^,  uq  does  not  have  a theoretical  background.  It 
is  derived  here  in  the  same  spirit  as  Cole9*  {Ref.  3)  law^of-the-wake. 

The  outer  variable,  however,  is  taken  as  y/0  rather  than  y/6  and  the  end 
result  is  entirely  different  from  the  law-of-the-wake.  Note  from  Eq.  (2) 

that  as  y+  ■+■  u!*  ■+■  tt/0.18.  ' Therefore,  according  to  Eq.  (1),  for  u+  to 

^ 1/2  + 

have  the  proper  limiting  value  of  (2/c^)  as  y + ®,  uq  must  behave  as 


as  y + ®.  Furthermore,  because  u.  gives  the  desired  result  for  small  y, 

+ + 1 
u must  behave  as  u + 0 as  y + 0.  Therefore,  the  form 

n r\  y ? 


] tan  1 

0.09 


(0.09y+) 


(6) 


is  considered.  Actually  g(y/9)  will  depend  on  cc,  H,  and  Re„,  but 
regardless  of  the  values  of  these  parameters,  the  function  g(y / 0 ) must 
behave  as  g(0)  = 0 and  g(™)  1. 


The  criterion  for  choosing  the  function  g(y/8)  was  that  it  be  a 

relatively  simple  function  that  does  a reasonable  job  of  describing  the 

trend  of  experimental  velocity  data.  The  experimental  data  trends  were 

investigated  by  solving  for  g(y/0)  from  Eq.  (6),  using  experimental  data 

+ + 

for  u , y , and  cf,  and  plotting  g(y/6)  versus  y/9.  An  example  of  the 
trends  is  given  In  Fig.  4 for  favorable  and  highly  adverse  pressure 


10 


AEDC -TR-77-79 


gradient  flows.  (These  data,  as  are  many  of  the  data  used  in  this 
report,  were  taken  from  the  tabulations  given  in  the  Stanford  Conference 
Proceedings  (Ref.  6).)  The  analytical  function  used  here  to  describe 
the  distribution  of  the  data  is 


(7) 


where  a and  b are  parameters  that  are  constant  for  a given  boundary- 

layer  profile  and  are  functions  of  c_,  H,  and  Re  . The  complete  velocity 

r tJ 

distribution  is  therefore  given  by 


where 


— L_  tan  1 
0.09 


(o.oV)  + "(fj /2 


tanh 


1/2 


-1  _ , and  JL_  „ ill 


(8) 


The  parameters  a and  b are  determined  by  satisfying  requirement 
number  (3) , which  is  that  velocity  profiles  similar  to  those  correlated 
by  von  Doenhoff  and  Tetervin  (Ref.  8)  be  recovered  away  from  the  wall  in 
the  outer  variable  y/6.  The  precise  profiles  established  in  Ref.  8 are 
not  recovered  because  of  the  fact  that  presumably  more  accurate  data 
were  obtained  subsequent  to  the  publication  date  of  Ref.  8 (1943),  and 
also  because  a Reynolds  number  effect  on  the  velocity  profiles  was 
reported  (see,  e.g. , Ref.  13)  subsequent  to  the  appearance  of  Ref.  8. 
Nevertheless,  the  idea  of  von  Doenhoff  and  Tetervin  is  what  is  important 
in  this  work,  and  that  is  that  turbulent  boundary- layer  velocity  dis- 
tributions can  be  correlated  (at  least  approximately  in  view  of  the 
Reynolds  number  effect)  in  a one-parameter  family  of  curves  of  u/u 

e 

versus  y/6,  with  each  curve  corresponding  to  a constant  shape  factor  H. 
Specifically,  a and  b are  determined  by  matching  the  distribution  u/ufi 
at  two  points,  y/0  - 2 and  y/6  » 5.  The  point  y/6  = 2 was  selected 
because  it  lies  relatively  close  to  the  wall  and  yet  is  outside  the 
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region  where  the  inner  solution  is  valid  (except  for  extremely  low 
Reynolds  number  flows) , and  y/0  = 5 was  selected  because  it  represents 
the  outer  region  and  is  a point  in  the  boundary  layer  where  the  velocity 
ratio,  u/ug,  is  relatively  constant  for  all  H as  demonstrated  in  Ref.  8. 


The  variation  of  u/ug(2)  (u/ue  at  y/0  = 2)  with  H is  illustrated  in 
Fig.  5a  by  the  experimental  data  tabulated  in  Ref. '6.  An  expression 
which  represents  these  data  reasonably  well  is  u/u  = 1.723  exp  (-0.6H) 
as  shown  in  Fig.  5a.  The  experimental  data  of  Ref.  6 are  for  the  most 
part  high  Reynolds  number  data,  and  data  for  Re0  >10^  were  used  in  Fig. 
5a.  The  effect  of  Reynolds  number  is  to  increase  the  ratio  u/u  (2)  as 
Re0  decreases  below  about  10  . Charts  were  established  by  Thompson 
(Ref.  13)  which  illustrated  this  trend,  and  these  were  used  as  a guide 
in  estimating  the  effect  of  Re0  on  u/ue(2).  An  empirical  expression 
which  accounts  for  the  variation  of  u/ue(2)  with  H and  Re0  is  given  by 


The  variation  of  u/u  (5)  (u/u  at  y/6  = 5)  with  H is  illustrated  in 
Fig.  5b,  again  using  the  experimental  data  of  Ref.  6.  Thompson’s  charts 
(Ref.  13)  indicate  a Reynolds  number  effect  for  these  u/ue(5)  data  also. 
However,  the  trend  suggested  by  the  charts  in  Ref.  13  does  not  appear  to 
be  entirely  substantiated  by  the  experimental  data  of  Ref.  6.  There- 
fore, the  Reynolds  number  effect  is  neglected  for  these  velocity  ratio 
correlations  at  y/0  = 5,  and  the  expression  used  to  correlate  the  data 
in  Fig.  5b  (which  includes  data  for  Re0  > 3,000)  is  given  by 


11 

u 

e 


(5)  = 0.87  + 0.08e 


-2.6UI-1.95)2 


(10) 


Note  in  Fig.  5a  the  data  extend  to  larger  H than  in  Fig.  5b.  The  data 
in  Fig.  5a  for  H > 2.9  are  measurements  made  by  Fraser  (Ref.  6)  in  a 
round  diffuser.  The  significance  of  this  is  that  these  axisymmetric 
internal  flow  measurements  correlate  well  with  the  two-dimensional  data 
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for  y/0  =°  2,  but  these  data  do  not  correlate  with  the  two-dimensional 
data  for  y/0  = 5.  Therefore,  for  axisyrametric  flows,  the  correlation 
given  by  Eq.  (10)  should  be  modified. 


Using  Eqs.  (9)  and  (10),  the  parameters  a and  b can  be  calculated 
as  described  in  detail  in  Table  1 . Table  1 is  simply  a summary  of  the 
solution  of  Eq.  (8)  for  a and  b at  the  two  points  y/0  = 2 and  5.  The 
curves  given  in  Fig.  .4  were  obtained  using  Eq.  (7)  with  a and  b calcu- 
lated in  this  manner.  The  distributions  given  by  Eq.  (7)  in  Fig.  4 are 
fair  representations  of  the  experimental  data,  particularly  in  view  of 
the  fact  that  any  discrepancies  in  Fig.  4 are  amplified  because  the 
experimental  data  were  referenced  to  u*  and  normalized  by 
i. e. , the  ordinate  in  Fig.  4 is  given  by  using  the  experimental  data  to 
evaluate 


/2_Y/2_  ir 
Vc£/  0.1 8J 


1 


Comparisons  of  Eq.  (8)  with  experimental  data  in  conventional  plots  of 

+ + 

u versus  y and  u/u£  versus  y/0  are  given  in  the  following  section. 


3.0  COMPARISONS  WITH  EXPERIMENTAL  DATA 


The  experimental  data  selected  for  comparisons  with  Eq.  (8)  are 
intended  to  represent  examples  weighted  to  flows  that  are  considered 
difficult  to  calculate  or  correlate,  as  opposed  to  flows  that  are 
readily  calculated  or  correlated.  For  example,  emphasis  is  given  to 
comparisons  with  boundary  layers  near  separation,  reattached  boundary 
layers,  and  nonequilibrium  boundary  layers.  In  addition,  compressible 
turbulent  boundary-layer  data  are  compared  with  Eq.  (8). 

The  much  used  experimental  data  of  Klebanoff  (Ref.  14)  and  Ludwieg 
and  Tillmann  (Ref.  6)  are  compared  with  Eq.  (8)  in  Fig.  6 in  the  vari- 
ables u and  y+.  In  Fig.  6 the  comparisons  involve  zero  and  mild  ad- 
verse pressure  gradient  flows  as  indicated  by  the  values  of  8 (where 
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0 = — -r*-  ) of  0.0  and  5.499.  Also,  a comparison  is  made  with  a pro- 
tw  ax 

file  (0  = 2.996)  which  comes  from  the  group  of  data  Ludwieg  and  Tillmann 
labeled  strong  adverse  pressure  gradient.  Actually,  this  0 = 2.996 
profile  was  measured  just  downstream  of  the  location  of  the  strongest 
pressure  gradient,  and  the  rate  of  change  of  pressure  gradient  was 
large.  In  all  three  examples  given  in  Fig.  6,  the  agreement  between  Eq. 
(8)  and  the  experimental  data  is  considered  good. 

In  Fig.  7 comparisons  in  u*"-y+  coordinates  are  made  between  Eq.  (8) 
and  two  adverse  pressure  gradient  flows  and  one  favorable  pressure 
gradient  flow.  The  data  of  Perry  (Ref.  6)  were  taken  in  a decreasing 
adverse  pressure  gradient  flow  and  resemble  (Ref.  6)  the  flow  in  a 
straight-walled  diffuser.  Stratford's  data  (Ref.  6)  were  taken  down- 
stream of  an  abrupt  onset  of  severe  positive  pressure  gradient.  Both 
Perry  and  Stratford's  data  are  out  of  equilibrium  (equilibrium  requires 
0 to  be  invariant  in  the  streamwise  direction)  and  Stratford's  data  are 
near  separation.  Bauer's  measurements  (Ref.  6)  were  made  in  water 
falling  down  a plate-glass  surface,  and  this  boundary  layer  was  near 
equilibrium.  The  agreement  between  Eq.  (8)  and  the  data  in  Fig.  7 is 
considered  good.  Note  that  there  is  little1,  if  any,  logarithmic  region 
remaining  in  Stratford's  profile. 

A difficult  test  for  Eq.  (8)  was  thought  to  be  that  provided  by  the 
experiments  of  Tillmann  (Ref.  6)  concerning  reattachment  and  recovery 
from  separation  at  a ledge  spoiler.  These  experiments  were  stated  in 
Ref.  6 to  be  among  the  best  available  concerning  this  problem.  Compari- 
sons are  made  in  Fig.  8 between  Eq.  (8)  and  four  of  Tillman's  profiles 
in  u/u^-y/6  coordinates.  The  two  profiles  nearest  the  ledge,  as  indi- 
cated by  the  circle  and  square  symbols  in  the  insert  in  Fig.  8,  are  also 
the  nearest  ones  to  the  ledge  which  were  reported  in  Ref.  6.  The 
agreement  in  Fig.  8 is  considered  reasonably  good. 


14 


AEDC -TR-77-79 


Of  practical  interest  for  aerodynamic  testing  are  the  measurements 
by  Newman  (Ref.  6)  of  an  airfoil  boundary  layer  proceeding  toward  separa- 
tion. Two  profiles  of  Newman  (Ref.  6)  are  compared  with  Eq.  (8)  in  Fig. 

9 in  u/u  -y/0  coordinates.  The  profile  of  the  last  axial  station  re- 
e 

ported  in  Ref.  6 is  used  in  Fig.  9 (8  = 182.776)  because  it  was  near 
separation,  and  the  other  profile  was  chosen  arbitrarily  and  has  a 8 of 
about  one  order  of  magnitude  less  (6  = 14.541).  The  agreement  between 
Eq.  (8)  and  the  experimental  profile  for  8 = 182.776  is  not  as  good  as 
the  agreement  between  Eq.  (8)  and  the  experimental  profile  for  6 “ 14.541. 
The  former,  however,  is  considered  acceptable  and  the  latter  extremely 
good.  Winter  (Ref.  15)  observed  that  a logarithmic  region  does  not 
exist  for  the  highly  adverse  pressure  gradient  profile  (6  = 182.776) 
shown  in  Fig.  9.  From  Fig.  10,  where  these  same  data  are  presented  in 
u+-y+  coordinates,  one  can  easily  see  Winter's  point.  The  analytical 
logarithmic  relation,  given  in  Fig.  10  for  reference,  contains  the 
constants  suggested  by  Coles  (Ref.  6).  Acceptable  agreement  is,  there- 
fore, obtained  between  Eq.  (8)  and  an  experimental  turbulent  boundary- 
layer  profile  that  has  no  logarithmic  region. 

Some  recent  measurements  by  Parikh,  Kays,  and  Moffat  (Ref.  16)  in 
low  Reynolds  number,  nonequilibrium,  adverse  pressure  gradient  boundary 
layers  are  compared  with  Eq,  (8)  in  Fig.  11.  The  profile  with  H = 2.002 
was  chosen  for  comparison  because  it  had  both  the  largest  value  of  H and  8 
of  all  the  profiles  reported  in  Ref.  16.  The  remaining  profile  in  Fig. 

11  was  chosen  to  represent  an  extremely  low  Reynolds  number  test.  The 
agreement  in  Fig.  11  between  Eq.  (8)  and  the  experimental  data  is  con- 
sidered reasonably  good.  However,  depending  on  c^  and  H,  the  data  in 
Fig.  11  for  Re^  ■ 1,220  are  near  the  conditions  where  the  point  y/B  = 2 
may  not  be  outside  the  region  where  the  inner  solution  is  valid.  In 
this  case  a different  point,  located  outside  y/9  = 2,  should  be  used  in 
place  of  the  point  at  y/0  = 2 for  determining  parameters  a and  b. 
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Compressible  turbulent  flows  are  now  investigated  from  the  point  of 
view  of  applying  Eq.  (8)  directly,  using  equivalent  incompressible 
values  for  the  variables  in  Eq.  (8),  and  assessing  the  quality  of  agree- 
ment with  experimental  data.  To  this  end,  the  transformation  theory  of 
Lewis,  Kubota,  and  Webb  (Ref.  17)  is  considered.  Lewis,  Kubota,  and 
Webb  suggested  the  coordinates 


and 


(11) 


(12) 


The  superscript  bar  ( ) notation  was  used  in  Ref.  17  and  refers  to  low- 
speed  flow.  Therefore,  y/6  is  taken  as  the  outer  variable  y/0  in  Eq. 

(8).  The  variables  defined  by  Eqs.  (11)  and  (12)  were  used  in  Fig.  2 of 
Ref.  17  to  correlate  a large  amount  of  high  Mach  number  constant  pressure 
turbulent  boundary- layer  data.  Figure  2 of  Ref.  17  is  used  as  Fig.  12 
of  this  report,  with  an  additional  curve  representing  Eq.  (8).  Equiva- 
lent incompressible  variables  for  Eq.  (8),  which  correspond  to  the 
compressible  data  in  Fig.  12,  were  determined  as  follows.  The  "law  of 
corresponding  stations"  (Ref,  20)  is  the  relation 


9 ? 

c.  Re„  = _i_l  (c.)  (ReJ 


f 6 p ii 


f c 


Q' 


(13) 


Assuming  a constant  pressure  normal  to  the  surface  and  a linear  varia- 
tion of  viscosity  with  temperature  as  was  used  in  Ref.  17,  Eq.  (13) 
reduces  to 


Cf  % “ (cfJe  (RVc 


(U) 


All  the  data  in  Fig.  12  are  for  5 < (cf)c(Re0)c  < 9.  Using  Eq.  (14) 

and  a mean  value  of  (cf)c(ReQ)c  - 7,  one  has  that  cf  ReQ  = 7.  Using  the 

skin  friction  correlation  of  Winter  and  Gaudet  (Ref.  21)  for  constant 

pressure  turbulent  boundary  layers,  and  the  mean  value  of  c-  Re„  = 7, 

£ 0 
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one.  has  that  = 0.0035  and  Re0  = 2,000,  approximately.  A representa- 
tive, incompressible  shape  factor  for  the  constant  pressure  turbulent 
boundary- layer  data  in  Fig.  12  is  taken  as  1.3  based  on  those  tabulated 
in  Ref.  21  for  compressible  flows  at  higher  Reynolds  numbers.  There- 
fore, Eq.  (8)  was  evaluated  for  the  equivalent  incompressible  values  of 
Cj  = 0.0035,  Re^  = 2,000,  and  H = 1.3.  The  agreement  in  Fig.  12  is 
considered  good,  particularly  so,  in  view  of  the  fact  that  hypersonic 
data  up  to  Mach  number  8.18  are  included.  Scatter  in  the  experimental 
data  in  Fig.  12  for  small  y/0  is  attributed  in  Ref.  17  to  the  lack  of 
accounting  for  viscous  dissipation  effects  in  the  transformation  theory. 
The  curves  marked  equivalent  incompressible  profiles  in  Fig.  12  were 
placed  there  by  Lewis,  Kubota,  and  Webb  (Ref.  17)  (recall  that  Fig.  12 
was  taken  from  Ref.  17).  Obviously,  Eq.  (8)  provides  a marked  improve- 
ment of  correlating  these  data  over  that  provided  by  the  equivalent 
incompressible  profiles  included  in  Fig.  12. 

Recent  solid  tunnel  wall  measurements  were  made*  in  the  Propulsion 
Wind  Tunnel  Facility  (PWT)  Acoustic  Research  Tunnel  (ART)  at  the  Arnold 
Engineering  Development  Center  (AEDC).  Two  of  the  measured  profiles 
representing  the  low  and  high  Mach  numbers  investigated,  Mg  = 0.20  and 
0.65,  are  compared  with  Eq.  (8)  in  Fig.  13.  Although  the  larger  Mach 
number  is  only  0.65,  it  is  sufficiently  large  that  compressibility 
should  be  taken  into  account,  particularly  with  regard  to  the  shape 
factor  H.  Compressibility  was  taken  into  account  by  first  reducing 
these  measurements  (which  were  pitot  pressure  measurements)  by  the  data 
reduction  method  described  in  Ref.  7 for  the  appropriate  free-stream 
conditions.  The  resulting  velocity  profiles  were  then  used  to  obtain 
the  corresponding  incompressible  values  of  H,  Re0,  and  cf  for  use  in  Eq. 
(8).  Relatively  good  agreement  between  Eq.  (8)  and  the  experimental 
data  is  demonstrated  in  Fig.  13. 


* These  measurements  were  made  by  Dr.  J.  A.  Benek,  ARO,  Inc. 
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Bound ary- layer  measurements  have  recently  been  reported  by  Altstatt 
(Ref.  22)  of  transonic  flow  over  a one-inch-high,  twelve- inch- long 
circular  arc  bump  placed  in  the  floor  (a  solid  floor  was  used)  of  the 
PWT  Aerodynamic  Wind  Tunnel  (IT)  at  AEDC.  A schematic  of  the  flow  and 
identification  of  the  location  of  the  measurements  are  given  in  the 
insert  in  Fig.  14.  The  data  Indicated  by  the  square  symbols  were  taken 
in  the  cusp  region  where  the  wall  is  faired  into  the  bump.  These  parti- 
cular measurements  were  made  by  Cline  (Ref.  22)  using  a Laser  Doppler 
Velocimeter  (LV)  system.  Compressibility  was  accounted  for  by  the 
method  described  above  for  the  ART  data,  except  that  velocity  data  were 
available  directly  from  the  measurements  because  an  LV  system  was  used. 
The  agreement  between  Eq.  (8)  and  the  experimental  data  in  Fig.  14  Is 
reasonably  good,  particularly  with  regard  to  the  measurements  in  the 
cusp  region. 


4.0  SUMMARY  OF  RESULTS 

An  analytical  expression  for  the  velocity  distribution  in  a tur- 
bulent boundary  layer  was  derived  and  shown  to  he  in  good  agreement  with 
the  experimental  data  considered  over  the  entire  domain  0 <_  y < «.  The 
analytical  expression  recovers  experimental  velocity  distributions  in 
the  region  near  the  wall,  matches  correlated  velocity  distributions  at 
y/0  =>  2 and  5,  and  gives  the  proper  limiting  velocity  as  y -*■  “.  The 
resulting  expression  gives  velocity  explicitly  as  a function  of  y and 
depends  on  properties  that  are  explicitly  defined. 

A distinguishing  feature  of  the  Inner  solution  results  is  that 
besides  the  good  agreement  between  the  inner  solution  and  velocity 
distribution  data  in  the  region  near  the  wall,  good  agreement  was  also 
obtained  with  Reynolds  stress,  turbulence  production,  and  turbulence 
dissipation  data.  To  the  extent  of  obtaining  this  agreement  with  experi- 
mental data,  the  inner  solution  is  therefore  consistent  in  that  the  same 
model  represents  all  these  mean  turbulence  quantities  simultaneously. 

An  additional  feature  is  that  the  expressions  for  these  turbulence 
quantities  are  exceedingly  simple. 
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The  single  major  difference  between  the  present  analytical  expres- 
sion of  the  turbulent  boundary  layer  velocity  distribution  and  preceding 
expressions  is  the  absence  of  an  explicit  logarithmic  term.  This  ab- 
sence of  a logarithmic  term  is  completely  contrary  to  years  of  tradition 
in  turbulent  boundary- layer  research,  but  was  the  primary  result  which 
eventually  led  to  the  development  of  a single  expression  that  provided  a 
reasonably  good  description  of  the  entire  (0  <_  y < «°)  velocity  profile. 
The  present  result  is,  no  doubt,  not  the  final  word  concerning  velocity 
distribution  expressions.  However,  the  important  result  of  the  present 
Investigation  might  be  the  demonstration  that  analytical  forms  other 
than  a logarithmic  distribution  (which,  incidently.  Coles  (Ref.  3) 
contends  is  empirical  because,  among  other  things,  of  the  requirement 
that  the  total  shear  stress  must  remain  constant  if  the  logarithmic  term 
is  to  be  derived  analytically,  and  this  requirement  contradicts  experi- 
mental observations  in  that  the  logarithmic  distribution  can  be  a good 
representation  of  the  actual  velocity  distribution  in  regions  where 
experiments  demonstrate  that  the  total  shear  stress  is  not  constant) 
should  receive  further  attention. 
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Figure  2.  Reynolds  stress  distributions  according  to  Eq.  (3)  and 
the  experimental  data  reported  by  Schubauer  (Ref.  12). 
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Figure  3.  Turbulence  production  and  dissipation  distributions 
according  to  Eqs.  (4)  and  {5}  and  the  experimental 
data  reported  by  Schubauer  {Ref.  12). 
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Figure  6.  Zero,  mild,  and  strong  adverse  pressure  gradient  boundary  layers  according 
to  Eq.  (8)  and  the  experimental  data  of  Klebanoff  (Ref.  14)  and  Ludwieg 
and  Tillman  (Ref.  6). 
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Figure  7.  Adverse  and  favorable  pressure  gradient  boundary  layers  according  to  Eq.  (8)  and  the 
experimental  data  of  Stratford  (Ref.  6),  Perry  (Ref.  6),  and  Bauer  (Ref.  6). 
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Figure  8.  Comparisons  of  the  experimental  date  of  Tillman 
(Ref.  6]  with  Eq.  {8}  for  boundary  layers  upstream 
and  downstream  of  separation. 
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Figure  9.  Airfoil  boundary-layer  flow  proceeding  toward 
separation  in  u/ue-y/0  coordinates. 
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Figure  10.  Airfoit  boundary-layer  flow  proceeding  toward  separation  in  u+-y+  coordinates. 
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Figure  11.  Low  Reynolds  number,  nonequilibrium,  adverse  pressure  gradient  turbulent  boundary  layers 
according  to  Eq.  (8)  and  the  experimental  data  of  Parikh,  Kays,  and  Moffat  (Ref.  16). 
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Figure  12.  Experimental  compressible  turbulent  boundary-layer 
data  in  transformed  coordinates  compared  to  Eq.  (8). 


35 


AEDC -TR-77-79 


Sjrm  cf  x IQ3  „ «eg  x IQ"3  ^ 

O 2.10  1.550  1.303  26.47  0.73 

C 1.37  1.673  1.502  40.96  0.61 

(cf,  H,  and  Re0  are  the  sane 

as  each  experiment.) 


Source 

Altstatt  (Ref.  22) 
Altstatt  (Ref.  22) 
Eq.  (8) 


Figure  14.  LV  data  from  the  Aerodynamic  Wind  Tunnel  (IT)  (Tunnel  IT) 
bump  experiments  (Ref.  22)  compared  to  Eq.  (8). 
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NOMENCLATURE 

a Parameter  in  Eqs.  (7)  and  (8) 

b Parameter  in  Eqs.  (7)  and  (6) 

Cp  Local  skin  friction  coefficient,  2t  / p 


8 

H 

M 

P 


u 


+ 

u 

+ 


X 


Function  defined  by  Eq.  (7) 
<5*/0 


Mach  number 


Static  pressure 


p u e/p 
e e e 


Mean  velocity  in  the  x direction 


<Wp..)1/2 


W V 


u/u 


Inner  solution  for  u 


Outer  expression  for  u+ 


Defined  by  Eq.  (11) 


Reynolds  stress 


Coordinate  along  body  surface 
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y 
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y 

y/e 

3 

6 

6* 

8 

y 

n 

p 

T 


Coordinate  normal  to  body  surface 

vA 

Defined  by  Eq.  (12) 

6^  d£ 

T dx 
w 

Boundary- layer  thickness 
Boundary-layer  displacement  thickness 
Boundary-layer  momentum  thickness 
Molecular  viscosity 
Coles'  profile  parameter 
Density 

Total  shear  stress 


SUBSCRIPTS 

c Compressible  value 

e Boundary-layer  edge  value 

w Wall  value 

aw  Adiabatic  wall  value 
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